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ABSTRACT 

Aims. We discuss whether the Gaussian is a reasonable approximation of the velocity distribution of stellar systems that are not spher- 
ically distributed. 

Methods. By using a non-Gaussian velocity distribution to describe the sources in the Large Magellanic Cloud (LMC), we reinvesti- 
gate the expected microlensing parameters of a lens population isotropically distributed either in the Milky Way halo or in the LMC 
(self lensing). We compare our estimates with the experimental results of the MACHO collaboration. 

Results. An interesting result that emerges from our analysis is that, moving from the Gaussian to the non-Gaussian case, we do not 
observe any change in the form of the distribution curves describing the rate of microlensing events for lenses in the Galactic halo. 
The corresponding expected timescales and number of expected events also do not vary. Conversely, with respect to the self-lensing 
case, we observe a moderate increase in the rate and number of expected events. We conclude that the error in the estimate of the 
most likely value for the MACHO mass and the Galactic halo fraction in form of MACHOs, calculated with a Gaussian velocity 
distribution for the LMC sources, is not higher than 2%. 
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1. Introduction 

Galaxies are complex, collisionless, gravitationally-bound sys- 
tems formed by secular gravitational self-interaction and col- 
lapse of its constituents. Significant progress has been made 
in both observational and theoretical studies developed to im- 
prove our understanding of the evolutionary history of galax- 
ies and the physical processes driving their evolution, leading 
to the Hubble sequence of galaxy type that we observe today. 
However, many aspects of their features, such as morphology, 
compositions, and kinematics, still remain unclear. In particular, 
it is not obvious how the velocities of their constituent (in par- 
ticular stellar) components can be described, because we cannot 
consider them to be isotropically distributed at any point. Little 
is known about the velocity distribution (VD) of the stellar pop- 
ulations of galactic components. While the distribution of stellar 
velocities in an elliptical galaxy is generally reasonably close to 
a Gaussian, analyses of the line-of-sight (l.o.s.) velocity distri- 
butions of disk galaxies have shown that these distribution are 
highly non-Gaussian ( |Binney & Merrifield 1998] >. 

Today, one of the most important problems regarding the 
composition of the Milky Way (MW) concerns the existence of 
dark, compact agglomerates of baryons in the Galactic halo, the 
so-called MACHOs (MAssive Compact Halo Objects). From the 
experimental point of view, several observational groups have at- 
tempted to detect these objects by performing microlensing sur- 
veys in the directions of the Large Magellanic Cloud (LMC), 
Small Magellanic Cloud, and M31. Two groups (MACHO and 
POINT-AGAPE) reported similar conclusions, despite the fact 
that they observed different targets (LMC and M31), that is 



roughly 20% of the halo mass must be in the form of MACHOs 
(Alcock et al. l2000l ICalchi Novati et al. 2005]) . However, the 
interpretation of their data is controversial because of the in- 
sufficient number of events detected, and the existing degen- 
eration among the parameters. Discordant results have been 
reported by other experimental teams (Tisse rand et al. 2 007; 
de Jong et al. 2004| i. 

Accurate theoretical estimates of the microlensing parame- 
ters, supported by statistical analysis, are fundamental to the in- 
terpretation of the experimental results. However, there are many 
uncertain assumptions in the adopted lens models. These uncer- 
tainties, that could lead towards an incorrect interpretation of the 
data, are mostly related to the shape of the individual galactic 
components and the kinematics of the lens population. 

One of the first problems to be raised by the scientific 
community concerned the shape of the Galactic dark halo. 
Unfortunately, information that can be extracted from observa- 
tions of high-velocity stars and satellite galaxies does not place 
strong constraints on its shape. In the absence of precise data, 
we are aided by computational models of the formation of galax- 
ies, which suggest that the dark halos are more or less spherical 
dNavarro et al. 1996i l. However, Griest dl991| ) showed that, re- 
ferring to MACHOs, the optical depth is relatively independent 
of assumptions about the core and cutoff radii of the MW halo. 
Sackett & Gould |(1993)| first investigated the role of the MW 
halo shape in Magellanic Cloud lensing, finding that the ratio 
of the optical depths towards the Small and Large Magellanic 
Clouds was an indicator of the flattening of the Galactic dark 
halo. Alcock et al. (2000) considered a wide family of halos, be- 
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sides the spherical one, ranging from a massive halo with a rising 
rotation curve to models with more massive disks and lighter ha- 
los. On the other hand, the problem related to the shape of the 
LMC halo was defined by Mancini et al. (2004). These authors 
also explored the consequences of different LMC disk/bar ge- 
ometries a part from the coplanar configuration. All these studies 
demonstrated that the estimate of the microlensing parameters 
were noticeably affected by the shape of the Galactic halo and 
the other Galactic components. 

For the kinematics of the lenses, the expression of their 
random-motion velocity was reanalyzed by Calchi Novati et al. 
(2006), who considered the LMC bulk motions including the 
drift velocity of the disk stars. This study indicates that the mean 
rotational velocity of the LMC stars is irrelevant to estimates of 
the MACHO microlensing parameters due to the preponderance 
of the bulk motion of the LMC. For the self-lensing, it is slightly 
significant for lenses located in the bar and sources in the disk. 
We emphasize that the VD of the LMC sources has always been 
modeled by a Gaussian. This assumption is just a first approxi- 
mation and was adopted for practical reasons. In this paper, we 
investigate whether this hypothesis is acceptable or not for differ- 
ent source/lens configurations or at least provide a quantitative 
measure of the effectiveness and accuracy of the Gaussian hy- 
pothesis. To achieve this purpose, we re-examine the framework 
of microlensing towards the LMC, and in particular we recalcu- 
late the number of expected events by assuming that the source 
velocities are no longer Gaussian distributed. Both the MACHO 
and the self-lensing cases are considered. 



2. Non-Gaussian velocity distributions 

If we consider a spherically symmetric distribution of stars with 
density p, then we can describe the dynamical state of the system 
by a distribution function of the following form 



F(E) 



„E/o- 2 



(27TO- 2 ) 



2^/2 



(1) 



where E — *P - v 2 /2 is the binding energy per unit mass, and *P 
is the relative gravitational potential (Binney & Tremaine 1987 ). 
It is well known that the structure of a collisionless system of 
stars, whose density in phase space is given by Eq. (Q]i, is iden- 
tical to the structure of an isothermal self-gravitating sphere 
of gas. Therefore, the velocity distribution at each point in the 
stellar-dynamical isothermal sphere is the Maxwellian distribu- 
tion /(v) = Ne^'- V ~l°" , which equals the equilibrium Maxwell- 
Boltzmann distribution given by the kinetic theory. If we con- 
sider a stellar system that is far from having a spherical distri- 
bution (for example, a galactic flattened disk, a triaxial bulge, 
or an elongated bar), we do not expect that it is correct to use 
a Maxwellian distribution to describe its velocity profile. In the 
same way, we must ask if it is correct or not to use a Gaussian 
shape f(v) ~ exp(-(v 2 /cr 2 )) to describe the l.o.s. or projected 
velocity profiles of non-spheroidal galactic components. We at- 
tempt to answer this question by using non-Gaussian VDs ob- 
tained by simulations. Numerical simulations of collapse and re- 
laxation processes of self-gravitating collisionless systems, sim- 
ilar to galaxies, are useful in identifying their general trends, 
such as density and anisotropy profiles. Two studies of the veloc- 
ity distribution function of these systems by means of numerical 
simulations were performed, showing that the velocity distribu- 
tion of the resultant quasi-stationary states generally becomes 
non-Gaussian (Iguchi et al. 120051 Hansen et al. 2006). 



2. 1 . Superposition of Gaussian distributions 

N-body simulations of different processes of galaxy formation 
were performed by Iguchi et al. (2005 ). As a result of their sim- 
ulations, these authors found stationary states characterized by a 
velocity distribution that is well described by an equally weighed 
superposition of Gaussian distributions of various temperatures, 
a so-called democratic temperature distribution (DT distribu- 
tion), that is 



/dt(v) = 4<! J-cre-^^M 



1 - Erf 



(jvn 



(2) 



where Erf(x) is the error function. The conclusion is that the DT 
velocity distribution is a universal property of self-gravitating 
structures that undergo violent, gravitational mixing. The origin 
of such universality remains, however, unclear. 

2.2. Universal velocity distribution 



Hansen et al. (2006) performed a set of simulations of controlled 
collision experiments of individual purely collisionless systems 
formed by self-gravitating particles. They considered structures 
initially isotropic as well as highly anisotropic. After a strong 
perturbation followed by a relaxation, the final structures were 
not at all spherical or isotropic. The VD extracted from the re- 
sults of the simulations was divided into radial and tangential 
parts. In this way, they found that the radial and tangential VDs 
are universal since they depend only on the radial or tangential 
dispersion and the local slope of the density; the density slope a 
is defined as the radial derivative of the density a = d In p/d In r. 
The points obtained by the simulations, which describe the uni- 
versal tangential VD, are described well by the following func- 
tional form, 
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where <T ta n is the tangential velocity dispersion, v tan is the two- 
dimensional velocity component projected on the plane orthog- 
onal to the l.o.s., while q and k are free parameters. Hansen et al. 
(2006) reported the universal tangential VD for three different 
values of the density slope a. Here we use the intermediate case 
where a equals -2. This VD has a characteristic break, where 
Vtan = l-6cr tan is taken to be the transition velocity. The low 
energy part is described by q — 5/3 and k = 0.93. By compar- 
ison, for the high energy tails, the parameters are q = 0.82 and 
k=l.3. 



3. Microlensing towards the LMC revisited 

Concerning the Hubble sequence type, the NASA Extragalactic 
Database considers the LMC as Irr/SB(s)m. The LMC is formed 
of a disk and a prominent bar at its center, suggesting that it 
may be considered as a small, barred, spiral galaxy. Different 
observational campaigns towards the LMC (MACHO, EROS, 
OGLE, MOA, SUPERMACHO) have been performed with the 
aim of detecting MACHOs. Among these, only the MACHO 
and EROS groups have published their results. The EROS col- 
laboration detected no events (Tisseran d et al. 20071 1. In contrast, 
the MACHO Project detected 16 microlensing events, and con- 
cluded that MACHOs represent a substantial part of the Galactic 
halo mass, but is not the dominant component (Alcock et al. 
2000). The maximum likelihood estimate of the mass m of the 
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lensing objects was w 0.5 M , whereas the fraction / of dark 
matter in the form of MACHOs in the Galactic halo was esti- 
mated to be ~ 20%. 

In the numerical estimates of the microlensing parameters, 
useful in studying the fraction of the Galactic halo in the form 
of MACHOs, a Gaussian shape velocity distribution is still com- 
monly used to describe the projected velocity distribution for 
the lenses as well as the source stars, although they are not 
spherically distributed (Jetz er et al. 20021 Mancini et al. 2004; 
lAssef et al. 20061 Calchi Novati et al. 2006 ). Here, our intention 
was to utilize the non-Gaussian velocity profiles described in the 
previous section for the sources, instead of the usual Gaussian 
shape, and show how the microlensing probabilities change ac- 
cordingly. As a concrete case, in Sect. 3.1 we analyzed two main 
parameters of the microlensing towards the LMC, the rate and 
the number of expected microlensing events generated by a lens 
population belonging to the MW halo as well as one belong- 
ing to the LMC itself. The results of our model were compared 
with the MACHO collaboration observational results (Alcock et 
al. [2000 ). Finally, the method of maximum likelihood is used 
in Sect. 3.2 to calculate the probability isocontours in the {m,f} 
plane. 



3. 1 . Microlensing rate and number of expected events 

We restricted our analysis by considering a homogeneous sub- 
set of 12 Paczyhski-like events taken from the original larger set 
B reported by MACHO; we did not consider the Galactic disk 
events, the binary event, and all candidates whose microlensing 
origin had been placed in doubt. In our calculations, we used 
the models presented in Mancini et al. (120041 ) to represent the 
various Galactic components: essentially an isothermal sphere 
for the Galactic halo, a sech 2 profile for the LMC disk, and a 
triaxial boxy-shape for the LMC bar. |van der Marel et al. (2 002) 
measured the velocity dispersion of the LMC source stars to be 
20.2 km/s. This measurement was completed as usual by a quan- 
titative analysis of the absorption lines in the LMC spectrum, by 
assuming a Gaussian form for the VD. In principle, to obtain an 
estimate of the velocity dispersion for a non-Gaussian distribu- 
tion, we have to repeat the same analysis of the LMC line profile 
by applying a non-Gaussian algorithm. To a first approximation, 
we ignored this subtlety and simply assumed that the dispersion 
of Gaussian and non-Gaussian VDs were equal. Fixing in this 
way the values of the velocity dispersions, we draw in Fig.[T]the 
universal tangential VD (dotted line) together with the TD dis- 
tribution (dashed curve), and a classical Gaussian profile (solid 
line). 

In general, the velocity of the lenses Ve consists of a global 
rotation plus a dispersive component. Since we assumed that the 
MW halo has a spherical form, we considered that the lenses 
are spherically distributed. In this case, the rotational component 
could be neglected, and at the same time we could safely con- 
sider the distribution of the dispersive component to be isotropic 
and Maxwellian ( |de Rujula et al. 1995[ ). This assumption was 
also supported by an analysis of the kinematics of nearly 2500 
Blue Horizontal-Branch Halo stars at |z| > 4 kpc, and with 
distances from the Galactic center up to ~ 60 kpc extracted 
from the Sloan Digital Sky Survey, where the observed distri- 
bution of l.o.s. velocities is well-fitted by a Gaussian distribution 
dXue et al. 20081 ). 

It is well known that the number of events N is the 
sum, N = 2 Afield i °f the number of events expected for 
each monitored field of the experiment defined to be A% e id = 
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Fig. 1. The velocity profiles used to describe the kinematics of 
the LMC stars: VDs derived from the simulations of Iguchi et 
al. (120051 ) (dashed curve) and Hansen et al. ( 120061 ) (dotted line). 
The solid black line represents a Gaussian profile. 
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Afield J (dr/d7 , E)£(rE)drE, where £fi e id is the field exposure, 
dF/dln is the differential rate with respect to the observed event 
duration, T E is the Einstein time, and S{T E ) is the detection ef- 
ficiency of the experiment. The differential rate is defined to be 
(Mancini et al. 2004; Calchi Novati et al. 12006] ) 
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where p s is the source density, /(v s ) represents the two- 
dimensional transverse velocity distribution of the sources, x 
is the ratio between the observer-lens distance D Q i and the 
observer-source distance D os , whereas p. is the lens mass in so- 
lar mass units. The normalization factor TV is the integral over 
the l.o.s. of the sources. The distribution dn(x)/dp represents 
the number of lenses with mass between p and /i + d/j at a 
given point in the Galactic component considered. Assuming 
the factorization hypothesis, we can write dn{x)/dp as the prod- 
uct of a distribution dno/dp depending only on p and the perti- 
nent density profile (de Rujula et al. 1995 Mancini et al. 2004) 
= m^^j where pe is the lens density. Concerning the 

functional form of dtio/dp, we supposed that for the lenses in the 
halo the mass function is peaked at a particular mass po, so that 
it could be described by a delta function dn^/dp = Sip - po)/po. 
For lenses in the LMC disk/bar, we utilized the exponential 
form d« /dju = Ap^e'^ 1 ^ dChabrier 20011) . where a = 3.3, 
po = 716.4, f3 = 0.25, whereas A is obtained from the normal- 
ization condition J Afx^e - ^^' dp = 1. 
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3.1 .1 . Lenses in the Galactic halo 

We calculated the differential rate of the microlensing events 
with respect to the Einstein time, along the lines pointing to- 
wards the events found by the MACHO collaboration in the 
LMC and for different values of fio- We used a Gaussian VD 
for /(v s ) as well as the non-Gaussian VDs, given by Eqs. (O and 
0). As yuo and the l.o.s. change, we did not observe any substan- 
tial reduction in the height of the distribution curve of the mi- 
crolensing event rate, and the corresponding expected timescale 
did not vary among the cases considered. With respect to the 
number of events, the situation did not change. Taking into ac- 
count the MACHO detection efficiency and the total exposure, 
we calculated the expected number of events, summed over all 
fields examined by the MACHO collaboration in the case of a 
halo consisting (100%) of MACHOs. Both in the Gaussian and 
the non-Gaussian case, we achieved the well-known result that 
the expected number of events is roughly 5 times higher than 
observed. 

3.1.2. Self-lensing 

We repeated the same analysis for the self-lensing configuration, 
that is where both the lenses and the sources are located in the 
disk/bar of the LMC. By varying the l.o.s., we found in general 
that the microlensing differential rate for the non-Gaussian case 
was higher than that for the Gaussian case. We noted that the ex- 
pected timescale also varied. Between the Gaussian and the non- 
Gaussian case, we also observed that the median value of the 
asymmetric distributions decreases of roughly 20%. Concerning 
the expected number of microlensing events, we estimated the 
same variation for sources with a non-Gaussian VD, that is an 
increase of roughly 20% from the value of ~ 1.2 events obtained 
with a Gaussian VD (Mancini et al. 120041 . 

3.2. MACHO Halo fraction and mass 

Following the methodology of Alcock et al. (120001) . namely the 
method of maximum likelihood, we estimated the halo fraction 
/ in the form of MACHOs and the most likely MACHO mass. 
The likelihood function is 




where N exp is the total number of expected events, and dF/dT^ 
is the sum of the differential rates of the lens populations 
(MACHOs, LMC halo, LMC disk+bar). The MACHO contri- 
bution is multiplied by /. The product applies to the N b s ob- 
served events. The resulting likelihood contours are shown in 
Fig- El where the estimate of the differential rate was performed 
using a Gaussian VD (solid line) and a universal VD (dashed 
line). The probabilities were computed using a Bayesian method 
with a prior uniform in / and m. A spherical isothermal distribu- 
tion was used to describe the lens density in the MW and LMC 
haloes. We found that the most probable mass for the Gaussian 
case is m = 0.60*033 wnere tne errors are 68% confi- 
dence intervals, and / = 23% with a 95% confidence interval of 
10% - 47%. We note that these values are slightly higher than, 
although fully compatible with, the original result reported by 
Alcock et al. (2000). The mismatch is due to some differences in 
the modelling and the fact that the set of the events considered is 
smaller. If we consider that the velocities of the stars in the LMC 
are non-Gaussian-distributed, the likelihood contours have mini- 
mal differences from those of the previous case. We note that the 
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Fig. 2. Likelihood contours for MACHO mass m and Halo frac- 
tion / for a typical spherical Halo. The contours enclose region 
of 34%, 68%, 90%, and 95% probability. The cross shows the 
maximum-likelihood estimate. Two different curves are shown 
for each contour according to the velocity profile adopted for the 
sources: a Gaussian shape (solid line) and a universal VD (gray 
dotted line). 

most significant variation is in the estimate of the lens mass, but 
that this is not higher than 2% for the 95% probability contour 
line. 

4. Discussion and conclusion 

We have investigated the limits of the validity of the Gaussian 
approximation used to describe the kinematics of a source pop- 
ulation in a microlensing context. This hypothesis, due to its 
practicality, has always been adopted without any check of its 
plausibility. We have remedied this deficiency in confirmation by 
an exhaustive analysis. To describe the motion of a stellar pop- 
ulation with a non-spheroidal distribution as correctly as pos- 
sible, we utilize two VDs (Sect. 2.1, Sect. 2.2) extracted from 
numerical simulations of collisionsless systems formed by self- 
gravitating particles. These VDs are substantially different than 
for a Gaussian one. We considered stars in the disk and bar com- 
ponents of the LMC and investigated their potential to be sources 
of lensing by transient lenses. In this framework, we recalculated 
the main microlensing parameters, including the MACHO halo 
fraction and the most likely value for the lens mass. 

For a configuration in which the lenses and sources belong 
to the target galaxy, we detected an increase in the differential 
rate of microlensing events towards the LMC when we used a 
non-Gaussian VD to describe the motion of its stars instead of 
a Gaussian one. This increase is reflected in the estimate of the 
number of expected events, which is roughly 20% higher than 
the 1.2 events found for the Gaussian case. 

The prediction for a halo that consists entirely of MACHOs 
is a factor of ~ 5 above the observed rates. The situation does not 
change in a noticeable way if we consider a non-Gaussian VD, 
since we have found that the number of events expected is practi- 
cally equal to the previous case. The results remain valid for both 
the DT and the universal VD. The main difference between the 
velocity dispersion of the LMC stars and the MACHOs, practi- 
cally neutralizes any possible variation due to the different shape 
of the VD of the sources. The maximum-likelihood analysis pro- 
vides values for m and / that are quite similar for the Gaussian 
and the non-Gaussian case. We conclude that the error in the es- 
timate of the most probable value for the MACHO mass as well 
as for the Galactic halo fraction in the form of MACHOs, cal- 
culated with a Gaussian VD for the LMC sources, is roughly of 
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the order of 1 - 2%. This fact implies that, in the study of the 
MW halo composition by microlensing, a Gaussian profile is a 
reasonable approximation for the velocity distribution of a sys- 
tem of source stars, even if they are not spherically distributed. 
On the other hand, for self-lensing, the Gaussian does not pro- 
vide a good description of the kinematics of a non-spherically 
distributed stellar population, in a similar way to the disk or the 
bar of the LMC. To ensure accurate microlensing predictions, it 
is thus necessary to replace the Gaussian VD by a more phys- 
ically motivated one, which takes into account the real spatial 
distribution of the source stars. 
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